Anomalous percolation and quantum criticality in diluted rare-earth nickelates 
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A microscopic model for the diluted spin-mixed compounds (RzYi-^BaNiOs (R=magnetic rare- 
earth) is studied using Quantum Monte Carlo (QMC). The ordering temperature is shown to be a 
universal function of the impurity concentration x and the intrinsic Ni-chain correlation length. An 
effective model for the critical modes is derived. The possibility of a quantum critical point driven 
by the rare earth-concentration and the existence of a Griffiths phase in the high dilution limit is 
investigated. Several possible experimental approaches to verify the results are put forward. 



An increasing number of low-dimensional quantum 
magnets are being realized in materials artificially mod- 
ified at nanometer scales. The use of layered structures, 
dilution and/or combination of magnetic species with 
different magnetic moment establish complex geometries 
that enhance the role of quantum correlations. The re- 
sulting competing ground states, quantum critical points, 
and coexistence of classical and quantum features are 
of great current interest in condensed matter and inter- 
disciplinary physics. The family of mixed-spin quantum 
antifcrromagnets (RzYi-^BaNiOs (R = magnetic rare 
earth) provides a unique experimental basis to discuss 
these issues. The x = compound is an excellent re- 
alization of the 1-dimcnsional (ID) 5=1 Heisenberg 
model, thanks to a linear-chain arrangement of the mag- 
netic Ni 2+ ions. Neutron scattering experiments [l| con- 
firmed a spin liquid behavior, and revealed a Haldane 
gap of about 10 meV in the magnetic excitation spec- 
trum. Inter-chain interactions, negligible for x = 0, 
can be controlled by substituting the intercalated non- 
magnetic Y 3+ by magnetic rare earths. For x = 1 the 
system orders antiferromagnetically at T/v=10-50 K, de- 
pending on R [2| ■ The spin dynamics is characterized by 
a coexistence and frequency-separation of semi-classical 
(spin waves) and quantum (Haldane-gap modes) excita- 
tions 0, 0, El By now, this phenomenon has been exten- 
sively studied experimentally , analytically 0, Q , and 
numerically |a,l3- 

The focus of the present Letter, however, is the role of 
randomness and disorder for < x < 1. Previous studies 
were aimed at understanding the large-x limit. Due to a 
finite intrinsic correlation length in Haldane spin chains 
(£c/i ~ 7 lattice units), for x w 1 disorder is effectively 
averaged out. The problem is then merely a rescaled ver- 
sion of x = 1 case 0, ■ For x < 1 /£ c /j the mean distance 
between magnetic rare earth centers becomes compara- 
ble to £. To date, very little is known of this regime, 
beyond the fact that it should represent qualitatively 
new physics. Below we develop a numerical approach to 
this interesting problem. Our calculations predict an un- 
usual percolation mechanism for (RzYi-^BaNiOs and 
related systems. We show that the ordering temperature 



is a universal function of x^ c h- In addition, we derive an 
effective model for the critical modes, discussing whether 
the system is ordered down to zero concentration x or 
can can support a quantum critical point. 

A microscopic model for the (R^Yi-^BaNiOs was 
uroposed in Refs. and studied in detail for x = 1 

7]. The corresponding spin Hamiltonian includes a part 
that describes the individual Haldane spin chains and 
an additional inter-chain interaction term. The former 
contains no disorder and is independent of the doping 
level x. In contrast, the interaction term involves spin 
operators for the intercalated rare earth ions, and directly 
reflects the effect of dilution. In the present work we will 
consider two variations of this model: 

Hi = i^c-c + J J]] Si 2j Si+l,2j (1) 

H 2 = H c-c ~ S i,2j S i+l,2j + r ^2 S i$j ( 2 ) 

ij ij 

with iJ c _ c = J c 2J ^i,2j+lS^ 2j+1 (Sf <2 j + 5? 2j - +2 ) 
ij 

Both Hamiltonians have a common inter-chain interac- 
tion term H c — C . e is a random variable representing 
the stochastic occupation of the inter-chain sites with 
moments Si^j+i corresponding to the rare earths with 
probability density p(e) = xS(e — 1) + (1 — x)S(e). As 
argued in Refs. |J, |7J, Kramers rare earth ions in the 
(RzYi-^BaNiOs structure can be modeled by 5 = 1/2 
pseudo-spins with anisotropic (Ising-type) R-Ni exchange 
interactions J c . 

The two Hamiltonians differ in their description of the 
gaped spin chains. In model the chains are S=l 
Heisenberg systems, as in the real material. This model 
is thus well suited for calculations of realistic ordering 
temperatures and direct comparisons with experimental 
data. However, model J5J allows more flexibility for a 
conceptual study of the problem. It incorporates quan- 
tum Ising spin chains in a transverse field (ITF) with 
r > J (paramagnetic gaped regime) ||. The ITF model 
is in the same universality class as 5 = 1 Haldane spin 
chains, and shares with it many common features 



The ground state in both cases is a spin liquid with ex- 
ponentially decaying correlations with the same asymp- 
totic power-law corrections. The spectrum is gaped, and 
the low-lying excitations are long-lived magnon-like state 
with a dispersion minimum near the ID AF zone-center. 
Though the symmetries of ITF and Haldane spin chains 
are different (the full SU (2) symmetry of the Heisenberg 
model is reduced to Zi in the ITF), it matters little in 
our case. Indeed, in the both complete Hamiltonians 
and El SU(2) symmetry is explicitly broken by the Ising 
nature of J c . The advantage of using the ITF model for 
our purpose is that it allows to tune the intrinsic in-chain 
correlation length £ c h = log (j,) through the entire range 
between the free-spin limit £ c /i = to chain-criticality 
£ c /j — > oo. This, in turn, enables a study of the fine in- 
terplay between the doping level x and the correlation 
length, and the resulting unique percolation behavior. 

Using an efficient Quantum Monte Carlo (QMC) clus- 
ter algorithm[loj[rH. we have first computed the critical 
temperature for a large range of transverse fields and con- 
centrations in (J2J using a periodic dilution pattern. In 
fig. QJi.) we plot, for different concentrations, the criti- 
cal temperature T c as a function of £ c /i A sharp crossover 
takes places when the temperature is of the order of the 
gap of the independent chain T — A c /j = 2(T—J) (dashed 
line) considered separate entity inside the whole 

system. It is remarkable that a purely quantum one- 
dimensional length scale is controlling the classical, two- 
dimensional transition at finite temperature. The key to 
understand the percolation mechanism in this system is 
the universal dependence of the critical temperature in 
the dimensionless variable xt; c h, which is the ratio be- 
tween the correlation length in the chain £ c h and the 
average distance between Nd 3+ ions For moderate 
values of the transverse magnetic field different concen- 
trations curves collapse when x approaches the low dilu- 
tion limit (see fig. EJ).) This suggests that the effective 
degrees of freedom that become critical interact with an 
effective coupling that decreases with the distance be- 
tween the Nd 3+ magnetic moments and the correlation 
length in the chains. To clarify the nature of such degrees 
of freedom we computed the critical exponents using the 
conventional finite-size analysis of the order parameter 
(see figHJl. We found that the exponents are compatible 
with the classical 2D Ising model as well for periodic as 
for disordered dilution (note that the critical behavior 
of the randomly diluted 2D Ising model differs from the 
one of the pure 2D Ising model only by logarithmic cor- 
rections) . This suggests that the building blocks of the 
effective critical model are classical Ising spins. 

To express this ideas in quantitative form we have 
derived a low energy effective model for both hamil- 
tonians using a canonical transformation. The central 
idea underlying this effective model is that the local mo- 
ments generated by the presence of the R ion between 
two chains are coupled due to the stiffness of the finite 





FIG. 1: a.) Critical temperature T c as a function of the cor- 
relation length in the independent ITF chain ^ c ^,for different 
concentrations of Nd 3+ x. The dashed line indicates the value 
of the gap of the independent chain (fcs — fi — 1). This en- 
ergy scale separates the two completely different behaviors at 
T c = A c h showing that the ID features of the gaped chains 
survive even at the critical point, b.) Universal dependence 
of the critical temperature on the dimensionless variable x£ c h- 
At moderate transverse fields there is a sizable range where 
the T c for different dilutions collapse to the critical tempera- 
ture of the 2D anisotropic Ising model (solid line). 
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FIG. 2: Finite-size analysis of the order parameter in the 
vicinity of the critical point for J c — 0.1, V = 1.25, x = 0.5 
and T c — 0.13 in the randomly diluted system. Here the data 
collapse is achieved with /3 = | and v = 1, the exponents of 
the 2D Ising model. 



chains connecting them. We first study the spectrum 
of the local hamiltonian of the Ni-R-Ni three-spin block: 

hs = JcsiMj-i + s ti+i) + r (^i-i + s t 3+ i)- Thc 

^-component of the central spin s z representing the mag- 
netic moment in the R 3+ ion is a good quantum num- 
ber. That symmetry breaks the block spectrum in two 
orthogonal sectors. The ground state is two-fold degen- 
erated and each ground state belongs to distinct sectors. 
We label these states |ff) and |JJ.) that have opposite val- 
ues for the expectation value of the total z component 
of block spin M z = ±(s z + S) where S = + p 2 )~^ 
is the expectation value of the two outer spin magne- 
tizations and fi — J C /2T. There is always a finite gap 
A B = ((J c /4) 2 +T 2 ) 1 / 2 to thc first excited state (level 
crossing of local states is forbidden by the symmetry). 
The local gap does not introduce new restrictions in the 
validity of our arguments because Ab > A c h and we are 
interested only in the regime T < A c h- All the formulas 
are valid also for the S=l system taking J c — > 2J C and 
r -> 

We consider now the dominant couplings in the low 
dilution limit: those between two effective local moments 
on neighboring sites along the direction of the chains i, i+ 
1. For model © this coupling is exclusively of the Ising 
type (there are no spin-flip process along the chains). 
The leading contribution to this coupling is given by the 
energy necessary to create a domain wall between blocks 
% and i + 1 , 



J\\ 



1 



where h» 



'II - JiS^Sfaj^+S&Sfaj^). Using the 
explicit form of the block eigenstates we can proof that 
the effective coupling is Jy = JS 2 + 0(J 2 /A B ), i.e. 
the product of the energy to create a domain wall in 
each chain times the expectation value of the magne- 
tization in the two outer spins of the block. In the 
case of S = 1 Heisenberg chains we have off-diagonal 
terms = J/2(S+ + Sf j+1 Sr +w )+ h.c. 

However, the sf are still good quantum numbers and 
this prevents spin-fliping in the effective Hamiltonian. 
This is the strongest argument in favor of the low-energy 
equivalence of models |QJ and An essentially identi- 
cal argument can be applied in the high dilution limit 
to compute effective coupling between two distant lo- 
cal magnetic moments obtaining now Ju = 2S 2 p where 
p = E(sj = s,Sk = -s) - E(sj — s,s k — s) E(sj,Sk) is 
the energy of the chain connecting the two blocks with 
spin Sj in site j and Sk in site k. In other words, p is 
the energy difference between anti-parallel and parallel 
spins at the end of the chain and it is directly related 
to the static stiffness of the (finite) chain which is non- 
zero even in the paramagnetic case. In both chain mod- 
els the asymptotic behavior is p ~ po cxp(— \k — j\/£ c h)- 
We can test numerically the validity of this model in 
the periodic diluted system. The critical temperature 




FIG. 3: (a) Normalized critical temperature T c {x) /T c (x — 1) 
for different models considered in the text. The arrows show 
the value of the inverse correlation length for each model, (b) 
Generic zero-temperature phase diagram of the randomly di- 
luted system described in the text. The paramagnetic region 
PM is a Griffiths-McCoy region. Note that the phase bound- 
ary of the ordered phase approaches V e ff — at x — with 
an essential singularity. The value of r c corresponds to the 
QCP of the pure system (see text). 



for the 2D anisotropic Ising model is given by the ex- 
pression sinh/3 c Ja; sinh /3 C J y = 1 where J x — J C S and 
J y = 2J 2 -^e"f and (3 C = 1/T C . The results (solid 
line) are compared with the numerical simulation in Fig. 
^d. The agreement emerges at T c << A# where thermal 
population of excited states above the doublet is strongly 
suppressed and the block is described by an Ising vari- 
able. Thc effect of thc quantum fluctuations in the micro- 
scopic models is integrated out in the effective classical 
model and can be formally separated into two parts: the 
reduction of the effective magnetic moment of the Ising 
spins |fr),|lj-) and the exponential decay of the couplings 
in the direction parallel to the chains. The most promi- 
nent feature of the model in the low dilution limit is the 
linear behavior of T c as a function of x at fixed ^ c h- This 
property survives in the more realistic randomly diluted 
case for a large range of concentrations (see Fig. |3]i)- 

It is interesting to consider whether the absence of per- 
colation threshold is a robust feature against the presence 
of weak spin-flip processes between the central spin R and 
its neighbors in the perpendicular direction that belong 



4 



to the chains, hf Y = J^ v /2s+ J (5r._ 1 + S,r. +1 )+ h.c. . 
A finite value for J£ v breaks the conservation law for s s z 
and establishes quantum coherence between |-fl-) and 
In leading order, it has the same effect of a transverse 
magnetic field acting on the spin blocks. If Jxy « Aj, 
we can compute within second order perturbation theory 
the matrix element r e g = (JJ- \H xy \ ff) = J| y /4A L . This 
T e f f should not be confused with the original transversal 
field r in the chains. The effective model at T=0 is now 
a 2D ITF in that has a quantum critical point (QCP) for 
j = 3.044(2) [l^. I n the original model this translates 
into a QCP governed by the concentration and therefore 
a finite percolation threshold, x~ x = £c/ik>g(£^) where 
a is a numerical constant. Due to the fact that the ran- 
dom dilution affects only inter-chain sites a percolating 
cluster of magnetic moments exists with probability one 
for any finite value of x. But unlike the homogeneously 
diluted S — 1/2 2D Heisenberg systems 13], it can be 
can be disordered by pure quantum fluctuations. 

In Fig. Et> we show the generic zero temperature phase 
diagram that we expect for the randomly diluted system: 
Since the critical value for the effective matrix element 
F e ff vanishes exponentially for x — > 0, strongly diluted 
(i< 1) experimental systems are most probably in the 
paramagnetic phase (PM). Nevertheless, since the sys- 
tem is disordered and close to a quantum phase transi- 
tion at T e ff(x) 1 one expects the region r e // < T c to be 
a Griffiths-McCoy phase [pj. Il5l llq . in which various 
physical observables display a singular behavior and ex- 
perimental measurements should show strong sample to 
sample fluctuations as we actually observed in the nu- 
merical simulations. The usual phenomenological argu- 
ment leading to, for instance, a diverging susceptibility 
in a diluted transverse Isin g sy stem close to but away 
from a QCP is the following [lTlllfij : Connected clusters 
with N spins appear with an exponentially small prob- 
ability p N = exp(-XN) (p being the probability for an 
occupied site), but first order perturbation theory tells 
us that for small transverse field strengths the gap of 
this cluster is also exponentially small in the number of 
connected spins A ~ exp(— aN), leading to an exponen- 
tially large tunneling time. This implies, at zero tem- 
perature, a power law behavior for the autocorrelation 
function [(of (t)af (0)] av ~ t~ x / a , and thus an algebraic 
singularity of, for instance, the local zero-frequency sus- 
ceptibility x(w = Q) ~ t- 1+x /°. These Griffiths-McCoy 
singularities are generic for a randomly diluted trans- 
verse Ising system (actually generic for any kind of dis- 
order). In our case the dilution is such that no isolated 
finite clusters exist, nevertheless parts of the system are 
more strongly coupled than others simply because there 
are more occupied inter-chain sites. Generic disorder, 
like random ferromagnetic bonds, also leads to Griffiths- 
McCoy singularities in non-diluted transverse Ising sys- 
tems, as was demonstrated numerically |l0j . hence we 



also expect them to be present in the system we were 
considerng here. 

In conclusion we would like to stress the experimen- 
tal implications of our theoretical findings. Our central 
predictions, particularly the universal scaling of the or- 
dering temperature, should be straightforward to verify 
on (R^Yi-ir^BaNiOs powder samples by bulk magnetic 
and calorimetric measurements or by neutron diffrac- 
tion. Conversely, the intrinsic chain correlation length 
can be effectively measured through measuring the x- 
dependence of the ordering temperature. The finer ef- 
fects discussed above, including those associated with 
strong disorder effects and algebraic singularities, can be 
directly probed in more involved single-crystal neutron 
scattering experiments. 
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